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SAFETY STOCK DETERMINATION 

The invention relates to a method for assigning to a set of products a set of corresponding 
continuous densities of demand distributions. More particularly it relates to a method for 
determining safety stock levels for several products and for controlling the corresponding stock 
level. 

TECHNICAL FIELD AND BACKGROUND OF THE INVENTION 

Warehouses are used as stock supply of on-sale items because of uncertainties of demand or 
because some suppliers only deliver large lot sizes. Apart from figuring out what product has to 
be ordered when and at what quantity, one of the most important problems in warehousing is to 
set the size of the safety stock correctly. The safety stock is supposed to cover the uncertainty 
of demand during the inbound lead time. If the safety stock level is too high then more financial 
resources are bound to the inventory than necessary. If, however, the safety stock level is too 
low, the supplier will more often run into out-of-stock situations and thus will not be able to 
satisfy the demand and thereby fail to meet the desired serviceability. 

Most of the software tools used in supply chain management and warehousing determine the 
level of safety stock for each of the products under the assumption that the daily demands are 
normally distributed or by simply setting a fixed safety stock level. However, the problem with 
this assumption is that in many cases the demand is far from being normally distributed. In 
particular, this is true for certain kinds of warehouses like for spare parts. In the end this results 
in incorrectly calculated safety stock levels, i.e. either the achieved service levels end up being 
too high and therefore the inventory level could have been reduced or the achieved service 
levels end up being too low, thereby failing to meet the desired serviceability. In both cases 
additional costs are generated. 

From Kohonen, Teuvo: "Self-Organizing Maps", 3rd edition, Springer Ser. in Information 
Sciences 30, Berlin 2001, a clustering algorithm is known. Additional information thereto can 
be found in Kohonen, Teuvo, "The self-organizing map", Neurocomputing 21, No. 1-3, 1-6 
(1998). 
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In Kullback, Leibler: "On Information and Sufficiency", in Annals of Mathematical Statistics 
22 (1951), pp 79-76, the Kullback-Leibler distance is described. 

In Scarf, "The Optimality of (s,S) Policies in the Dynamic Inventory Problem", in Arrow, 
5 Suppes (Ed.), Mathematical Methods in the Social Science", Stanford University Press, 
Stanford 1960" the (s,S) replenishment method is described. More about that can be found in 
Zheng, Y.-S. "A simple Proof for the Optimality of (s,S) Policies in Infinite Horizon Inventory 
Systems", Journal of Applied Probability 28, 802-810 (1991) 

In Schneider, H., "Effect of Service-Levels on Order-Points or Order-Levels in Inventory 
10 Models", International Journal of Productions Research 19, 615-631 (1981) service level 
definitions are described. 

OBJECT AND ADVANTAGES OF THE INVENTION 

According to an aspect of the invention as set forth in the claims a method for assigning to a set 
15 of products a set of corresponding continuous densities of demand distribution is provided . 

The proposed method has the advantage of offering a safety stock level with higher precision, 
i.e. offering a product stock controller a means to get closer to his target to minimize the overall 
safety stock value under the condition of maintaining a given service level. The precision is 
particularly achieved utilizing a clustering step that on one hand reduces the number of 

20 determinations in a subsequent selection step and that on the other hand is able to introduce a 
cross-product statistical data enhancement effect, i.e. improve the demand statistics especially 
for sporadically demanded products. This allows to use more density models in the selection 
step and thereby achieve better precision in that step. Using the same number of density models 
for direct comparison with the products 1 densities would lead to an excessive uneconomic 

25 determination workload. 
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According to another aspect of the invention a method is provided that further encompasses a 
safety stock determination wherein with a given service level from the continuous demand 
densities a product safety stock level is determined. 

According to yet another aspect of the invention the method can be enhanced by a product 
5 order step wherein via an output device an order for a product is initiated towards a 
product-providing entity, e.g. via a communication medium like a communication network, in 
the event that a stock level of that product falls below a reorder point. Thereby an automated 
product stock control is achievable. 

The method can further comprise a clustering correction step wherein for each product a 
10 distance between its continuous density and its normalized discrete demand density, is 
determined and wherein for those of said products whose distance exceeds a predetermined 
threshold value, in a substitute-selection step for the corresponding normalized discrete demand 
density out of said predetermined set of continuous model densities a substitute continuous 
density is selected, which is defined for that product as its continuous density. That cluster 
15 correction step is advantageous because it is able to perform a distance check between the 
determined continuous densities and the normalized discrete demand densities and perform a 
correction should the distance be found too large. This improves the precision of the resulting 
assigned continuous densities and thereby reduces the error in the final safety stock levels. 

In the conversion step the demand time series can be substantially segregated from demand 
20 trend and periodic demand. Segregating the demand time series from demand trend and 
periodic demand is another way of increasing precision, since the parts that are removed would 
lead to errors in the resulting assigned continuous densities and hence in the final safety stock 
levels. 

In the normalization step normalization information can be maintained for the adjustment step. 
25 Maintaining the normalization information for the adjustment step facilitates the method since 
that information can be directly used there and need not be constructed in a different manner. 

If for each of the clusters an initial cluster representative is provided, and if the clustering step 
comprises for each of the normalized discrete demand densities a distance determination step 



CH9-2001-0064 

-4- 

wherein for each said initial cluster representative its distance from the normalized discrete 
demand density is determined, and a modification step wherein the one initial cluster 
representative is modified towards the normalized discrete demand density to which it has the 
smallest distance, said modified cluster representative being thereafter defined as new initial 
5 cluster representative. Using a difference estimation step and a modification step in the 
clustering step provides advantageous because this method is a quick and precise way of 
clustering. Furthermore the Kohonen algorithm can be utilized therefor. 

If the distance determination step and the modification step are repeated one or more times for 
each of the normalized discrete demand densities, starting with the initial cluster 
10 representatives that are present after the previous last modification step, this provides the 
advantage that the clustering result will be preciser which again leads to less errors in the 
resulting assigned continuous densities and hence in the final safety stock levels. 

Selecting the initial cluster representatives that are present after the last modification step in the 
clustering step as the cluster-representative discrete density brings the advantage that since the 
15 initial cluster representatives are already present, no additional steps for the selection have to be 
performed. This hence reduces the time the method takes. 

If in the clustering step a said cluster representative is substantially not modified when its 
distance to those normalized discrete demand densities that belong to the corresponding cluster 
is below a predetermined value, unnecessary representative movements are avoided, which on 
20 one hand reduces the clustering time and on the other hand avoids an exacerbation of the 
resulting clusters. Adding a condition that this cluster need have a number of thereto-belonging 
demand densities, also called members, that is bigger than a predetermined number, would 
introduce the advantage that only clusters of a predetermined minimum size are handled this 
way, again improving the clustering result. 

25 If the predetermined set of continuous model densities contains the model densities with its 
parameters in an undetermined form and the parameters for the clusters, are determined in the 
selection step, this brings in the advantage that the selection step becomes less 
storage-intensive since less model densities need be stored. 
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If in the normalization step the demand value for 0 pieces of the products in the discrete 
demand densities is rescaled inversely to the demand values for more than 0 pieces of the same 
product, a preciser clustering for products with sporadic demand, also referred to as 
slow-moving items, is the advantageous result. 

If in the selection step the demand value for 0 pieces of the products is suppressed, the 
selection result is less distorted by that demand value, which improves the resulting selection. 

If the method further comprises a safety stock determination step wherein with a given service 
level from the continuous demand density for one or more of the products a product safety 
stock level is determined, a stock controller gets preciser information about how many pieces 
of the product should be kept on stock in order to guarantee the service level. 

Performing a product order step wherein via an output device an order for a said product is 

initiated towards a product-providing entity, in the event that a current stock level, of that 

product falls below a minimum quantity that is the sum of the average of the discrete demand 

density for said product during the product's lead time, and the determined safety stock level for 

said product provides advantageous since it realizes an automated stock control for the 

products, thereby decreasing the probability of out-of-stock situations. 

« 

The above methods are implementable in form of program code means for performing the 
method. 

SUMMARY OF THE INVENTION 

The invention is directed to a method for assigning to a set of products a set of corresponding 
continuous demand densities. In a conversion step for each product its demand time series is 
converted into a discrete demand density. In a normalization step the discrete demand densities 
are transformed into normalized discrete demand densities. In a clustering step each of the 
normalized discrete demand densities is assigned to a cluster and for each cluster a 
cluster-representative discrete density is determined. In a selection step for each 
cluster-representative discrete density out of a predetermined set of continuous model densities 
a cluster-representative continuous density is selected. In a parameter-determination step for 
each product for its cluster-representative continuous density product-individual density 
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parameters are determined under use of which for each product a continuous density is 
determined. In an adjustment step for each product the continuous density is adjusted to form 
therefrom the continuous demand density. 

DESCRIPTION OF THE DRAWINGS 

5 Examples of the invention are depicted in the drawings and described in detail below by way of 
example. It is shown in 

Fig. 1 a block diagram, specifying the different steps and results of a density assignment 
process. 

Fig. 2 a schematic view of a computer in which from given product-specific input parameters a 
10 stock-alert signal or an automatic product order are issued. 

All the figures are for sake of clarity not shown in real dimensions, nor are the relations 
between the dimensions shown in a realistic scale. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

In the following, the various exemplary embodiments of the invention are described. 

15 For sake of simplicity, in the following, when referring to a "density", the density of the 
demand distribution is meant. When the index "i" is used herein for a specific reference 
indicator this is meant as a variable that could take any value of the specific values used 
elsewhere in this description for the same reference indicator. 

In figure 1 a block diagram, specifying the different steps and results of a density assignment 
20 process is shown. 

For each of a multitude of products P b P 2 , P 3 , P 4 , a corresponding demand time series H(P V ), 
H(P 2 ), H(P 3 ), H(P 4 ) is given. The demand time series can be a historical demand data, e.g. the 
demand development for the respective product Pi, P 2 , P3, P4 watched over a specific period of 
time. The time is herein quantified into time units, e.g. days. The demand is herein quantified 
25 into demand units, i.e. the smallest measured demand unit, to allow comparison of different 
product types. E.g. for screws a package of 1000 pieces can be regarded as smallest traded and 
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demanded unit, whereas for orange juice, 1 liter will be regarded as smallest traded and 
demanded unit. These demand units will in the following be referred to as demand quantities. 

In a conversion step 10 for each product Pi, P 2 , P 3 , P4 its demand time series H(Pi), H(P 2 ), 
H(P 3 ), H(P 4 ) is converted into a corresponding discrete demand density D P i d , D P2d , D P3d , Dp 4d by 
5 a binning process that defines a finite set of demand bins into which the actual demands are 
put. Such a bin could e.g. provide for a demand from 1 to 5 pieces belonging into a first bin, 6 
to 10 pieces belonging to a second bin, a.s.o. Each of these discrete demand densities D P i d , D P2d , 
D P3d , D P4d hence shows for its corresponding product Pi, P 2 , P 3 , P4, how often a specific amount 
of that product P,, P 2 , P 3 , P 4 is demanded. This discrete demand density D P i d , D P2d , D P3d , D P4d is 

10 free from demand trend and periodic demand. The purging can be done in the conversion step 
10, but can also be done before, i.e. the demand time series H(Pi), H(P 2 ), H(P 3 ), H(P 4 ) can also 
be already provided in a form that is free from demand trends and periodicities. "Free" herein 
means substantially free, since the purging in reality can not guarantee a 100% removal of trend 
and periodicity. Thus in the discrete demand density D P i d , D P2d , D P3d , D P4d substantially the 

15 stochastical component remains. Since demand trend and periodic demand are those demand 
components that can be more easily predicted, the safety stock need not cover it. The 
remaining, stochastical component is hence the one for which the safety stock is provided for 
the corresponding product Pi, P 2 , P 3 , P 4 . 

As next step follows a normalization step 20 wherein the discrete demand densities D P i d , D P2d , 
20 D P3d , D P4d are transformed into normalized discrete demand densities D N P i d , D N P2d , D N P3d , D N P4d . 
The normalization provides for a scaling to a norm value for the average and a norm value for 
the variance. The typical norm value for the average is 0 and the typical norm value for the 
variance is 1. The normalization to the norm value for the average effects a shifting by a 
shifting-factor s a , whereas the normalization to the norm value for the variance effects a scaling 
25 by a scaling-factor s v . These two factors s a , s v can be maintained for each normalized discrete 
demand density D N P i d , D N P2d , D N P3d , D N P4d> respectively the corresponding product Pi, P 2 , P 3 , P 4 , 
as adjustment factor set v ss i, v ss2 , v ss3 , v ss4 for a later adjustment step 70. The normalization has 
the effect that in the subsequent step, the densities D N P i d , D N P2d , D N P3d , D N P4d , that are compared 
have a common scale which effects a preciser relative positioning of the densities D N P i d , D N P2d , 
30 D N P3d , D N P4d , with respect to each other. More precisely, the distances dc(Xi, D N Pid ), as explained 
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below, will be more accurate which leads to a fairer clustering, i.e. those densities D N P , d , D N P2d , 
D N p 3d , D N P4d , which have the same shape will eventually be assigned to the same cluster. 

As next step follows a clustering step 30 wherein each of the normalized discrete demand 
densities D N P i d , D N P2d , D N P3d > D N P4d is assigned to a cluster Ci, C 2 and wherein for each cluster Ci, 
C 2 a cluster-representative is determined which is a cluster-representative discrete density Rfc. 
The number of clusters Ci, C 2 is predefined and can be determinable for this method. 

The clustering step 30 comprises for each of the normalized discrete demand densities D N P i d , 
D N P2d , D N P3d , D N P4d two steps, a distance determination step and a modification step. First for 
each of the clusters Ci, C 2 an initial cluster representative R ic , also referred to as initial cluster 
representative Xi, X 2 , is provided. The initial cluster representative Xi, X 2 need not be a 
normalized density since in the subsequent steps the quality of being normalized does not play 
a role. Normalization is here only of value for the relative positions of the discrete demand 
densities D N Pld , D N P2d , D N P3d , D N P4d . 

In the distance determination step, for all initial cluster representatives Xi, X 2 their distance 
d c (X i? D N Pid ) from one of the normalized discrete demand densities D N P i d , D N P2d , D N P3d , D N P4d is 
determined. The one normalized discrete demand density D N PJd , D N P2d , D N P3d , D N P4d can be 
selected according to any selection mechanism, be it a simple sequential scheme, or a random 
selection. 

As distance herein the Kullback-Leibler distance can be used, by defining each cluster 
representative X 1? X 2 and each normalized discrete demand density D N P i d , D N P2d , D N P3d , D N P4d as 
a point in an n-dimensional Kullback-Leibler space, wherein n is the number of quantity bins 
used in the conversion step 10. The position of the normalized discrete demand densities D N P i d , 
D N P2d , D N P3d , D N P4d in the Kullback-Leibler space depends on their average and variance, which 
due to the normalization step 20 has been unified, thereby leading to a common determination 
basis and a correct relative positioning within the Kullback-Leibler space. The positive effect of 
the normalization step 20 is present for any type of clustering involving distance-determination. 

The start position of the initial cluster representatives Xi, X 2 can be preselected, but it is found 
that performing a number of iterations, as is described below, the start position becomes 
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practically irrelevant. Therfore a start position of the initial cluster representatives Xi, X 2 can be 
selected, that is independent from the type, number or density of the products P b P 2 , P3, P4. The 
method can be designed to renounce a possibility to determine start positions and thereby 
require less user interaction by a stock controller. 

5 Next follows a modification step. In this step among the found distances dc(X i? D N P i d ) the 
smallest distance and the corresponding pair of initial cluster representative Xi and normalized 
discrete demand density D N Pid is determined. Then that initial cluster representative Xi is 
moved, i.e. modified, towards that normalized discrete demand density D N P i d , if the distance 
dc(Xi, D N Pid ) is not already below a predetermined maximum distance level, e.g. 0. The 
10 Kohonen clustering algorithm can be used to perform that function, wherein the cluster 
representatives Xi will be what in the Kohonen algorithm is referred to as centers. 

The distance determination step and the modification step are performed for each of the 
normalized discrete demand densities D N P i d , D N P2d , D N P3d , D N P4d . In the case of only a single run 
of the distance determination step and the modification step, it is useful to guarantee that every 
15 normalized discrete demand density D N P i d , D N P2d , D N P3d , D N P4d is once handled this way. In case 
of multiple iterations, a random selection is deemed to satisfy this criterion by effect of its 
statistical uniform distribution. In that case, a tracking of the use of each normalized discrete 
demand density D N P i d , D N P2d , D N P3d , D N P4d is not necessary. 

The clustering step 30 hence has the effect that the clusters Ci, C 2 and their representatives Xi, 
20 X 2 are defined and that furthermore it is defined which of the normalized discrete demand 
densities D N P i d , D N P2d , D N P3d , D N P4d belongs to which of the clusters Ci, C 2 . 

The clustering step 30 can be repeated in that in the distance determination step for each 
modified cluster representative Xi, X 2 its distance dc(X i5 D N Pid ) from the normalized discrete 
demand densities D N P]d , D N P2d , D N P3d , D N P4d is determined, and in the modification step each 
25 modified cluster representative X b X 2 is modified towards the normalized discrete demand 
density D N P i d , D N P2d , D N P3d , D N P4d to which its distance d c (Xi, D N Pid ) is smallest. This repetition, 
also called iteration, can be done several times, referred to as epoque in the Kohonen algorithm. 
Thereby the accuracy of the clustering step 30 is increased. This means that the final clusters 
Ci,C 2 are more precisely defined, meaning that the modified cluster representatives Xi, X 2 lie in 
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a closer distance to the normalized discrete demand density D N pid, D N p2d, D N P3d, D N P4d that 
belong to its cluster G, C 2 . 

It has been found that irrespective of the number of products P b P 2 , P 3 , P4 and even irregard of 
the density types thereof, a number of between 30 and 50 iterations are sufficient to achieve a 
5 sufficient stability. Generally it can hence be stated that it is advantageous to perform more 
than 30 iterations. Since with the number of iterations the execution time rises, it can also be 
stated that less than 50 iterations are advantageous. As stability criterion the average distance 
dc(Xj, D N Pid) between the modified cluster representatives X u X 2 and the normalized discrete 
demand density D N P i d , D N P2d , D N P3d , D N P4 d that belong to the respective cluster G, C 2 , can be 
10 used. Also an automatism can be used that stops performing further iterations if a certain 
stability of the sample locations has been reached, i.e. the average distance dc(Xi, D N Pid ) has 
fallen below a predetermined maximum distance value. 

The distance determination step can be executed under use of different distance determinators, 
thereby providing several values for the distances dc(Xi, D N P i d ). In such a case the modification 

15 step is executed using the smallest value of the distance dc(Xi, D N Pid ). The use of different 
distance determinators allows to avoid situations in which the distance d c (Xj, D N Pid ) that is 
calculated and determined does not correctly reflect the true relation between the modified 
cluster representative Xi, X 2 and the normalized discrete demand densities D N P i d , D N P2d , D N P3d , 
D N P4d . The different distance determinators typically perform different distance-determination 

20 algorithms that hence also result in different values for the distances d c (Xi, D N Pid ). It is 
assumable that the distance determinator with the smallest value of the distance d c (Xi, D N P j d ) is 
the most accurate and useful. The term "distance" is to be understood as a metric in the 
mathematical sense. 

Finally a cluster-representative discrete density Rfc is determined for each cluster G, C 2 . The 
25 simplest selection is to select the modified cluster representatives Xi, X 2 that exist after the last 
modification step as the cluster-representative discrete density R^. This selection has the 
advantage that no additional algorithm for selecting the cluster-representative discrete density 
Rfc has to be performed. It is also assumable that after running through the modification step, in 
particular after several iterations, the modified cluster representatives X b X 2 are the best 
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possible choice for the cluster-representative discrete density Rfc. Since the 
cluster-representative discrete density R rc is only used for a calculation step, it need not be any 
of the normalized discrete demand densities D N P i d , D N P2d , D N P3d , D N P4d . The clustering step 30 
has a very valuable advantage in that by the clustering products Pi, P 2 , P 3 , P4 are grouped that 
5 not necessarily have any correlation to each other. In particular this means that those products 
Pi, P 2 , P 3 , P4 which have only few demand occurrences in the period of time used for the 
demand time series H(Pi), H(P 2 ), H(P 3 ), H(P 4 ), in the following referred to as slow-moving 
items, profit from those products Pi, P 2 , P 3 , P4 that have a similar normalized discrete demand 
density D N P i d , D N P2d , D N P3d , D N P4d , but which have a higher number of demand occurrences in the 
10 same period of time. This increases the accuracy of the subsequent selection step 40 for the 
products Pi, P 2 , P 3 , P 4 with only few demand occurrences. 

The purpose of clustering the data is hence twofold. First, it enhances the quality of the 
succeeding steps by enlarging the amount of data that is being analyzed. This is especially true 
for slow-moving items because those have only little statistical data and the clustering effects 
15 the use of statistical data of all products Pi, P 2 , P 3 , P 4 in that cluster. Second, the subsequent 
steps need only be performed on the cluster's representatives X b X 2 . 

The clustering step 30 can be further improved by introducing an exceptional rule: In the 
clustering step 30, the cluster representative Xi, X 2 is substantially not modified when its 
distance d c (Xj, Pi) to the normalized discrete demand densities D N P i d , D N P2d , D N P3d , D N P4d that 

20 belong its cluster G, C 2 , is below a predetermined distance value. By giving the distance value 
for the distance dc(Xi, D N Pid ), being a maximum value, a cluster representative X b X 2 is 
substantially prevented from being moved away from these normalized discrete demand 
densities D N P i d , D N P2d , D N P3d , D N P4d , to which it is so close. This reflects the idea that in the case 
one of the cluster representatives X b X 2 is very close, i.e. the distance dc(Xi, D N Pid ) being below 

25 the predetermined distance value, to the corresponding normalized discrete demand densities 
D pid, D N P2d , D N P3d , D N P4d , then the cluster representative Xi is to be regarded as a good 
representative and hence shall not be modified, i.e. not by some effect of the clustering step 30 
be moved away from its good position. This cluster representative Xi, X 2 is practically frozen 
to its position, avoiding unwanted exacerbation of its position and even reducing the 

30 determination time in the modification step. Additionally the condition that the cluster d, C 2 
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should have a minimum size, i.e. a minimum number of thereto-belonging normalized discrete 
demand densities D N P i d , D N P2d , D N P3d , D N P4 d, may be added, to ensure that the cluster freezing is 
only performed for clusters Ci, C 2 of significant size. 

As next step follows a selection step 40 wherein for each cluster-representative discrete density 
5 R rc out of a predetermined set L D of continuous model densities D m Mc a cluster-representative 
continuous density D m Cc is selected. This is hence the step that goes from the discrete densities 
to the continuous densities. The task is to find among the continuous model densities D m Mc one 
that fits best the respective cluster-representative discrete density Rfc. The set L D should 
comprise model densities of different basic types like Gauss, Gamma, Poisson, Weibull, 
1 0 Logistic, and other parametrized density families. 

The predetermined set L D of continuous model densities D m Mc here contains the model densities 
D m Mc with their parameters p in an undetermined form. This has the advantage that the set need 
not contain the model densities D m Mc in a lot of possible parametrically determined forms which 
would unnecessarily increase the need for storage and determination time for going through all 

15 of them. Instead, in the selection step 40, the mere shape of the model densities D m Mc is 
compared to the cluster-representative discrete density Rfc to identify a best fit. An algorithm 
suitable for this operation would be the method of moments or Maximum Likelihood 
Estimators (MLE). In the course of the selection step 40 it is possible to have some of the 
parameters p of the model densities D m Mc already determined once they are selected for each 

20 cluster-representative discrete density Rfc. This is in particular useful for shape parameters. 

After the selection step 40 follows a parameter-determination step 50 wherein for each product 
Pi, P 2 , P 3 , P4 for its cluster-representative continuous density D m Cc product-individual density 
parameters p(Pi) are determined under use of which for each product Pi, P 2 , P 3 , P4 a continuous 
density D m P i c , D m P2c , D m P 3c, D m P4c is determined. This step hence leads back from the clusterwise 
25 handling to the productwise handling. 

The cluster-representative continuous density D m Cc that has been determined for each cluster Ci, 
C 2 is now for each product Pi, P 2 , P 3 , P4 that belongs to the respective cluster Ci, C 2 compared 
to the corresponding normalized discrete demand density D N P i d , D N P2d , D N P3d , D N P4 d and by 
means of that comparison the product-individual density parameters p(Pi) for that product Pi, 
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P 2 , P 3 , P 4 are determined. This leads to the continuous density D m Pt c, D m P2c , D m P3 c, D m P4 c for each 
product Pi, P 2 , P 3 , P4. It is possible to use the already determined parameters p, should those 
already have been determined in the selection step 40. 

At this point in the method the continuous densities D m P i c , D m P2c , D m P3c , D m P4c for each product 
5 Pi, P 2 , P3, P 4 have been determined. On the basis of those densities D m P i c , D m P2c , D m P3c , D m P4c the 
safety stock level $ for all the products Pi, P 2 , P 3 , P 4 can be determined. It is however possible 
that due to algorithmic uncertainties, the found continuous density D m P i c , D m P2c , D m P3c , D m P4c 
does not accurately enough reflect the normalized discrete demand density D N P i d , D N P2d , D N P3d , 
D N P4d . A typical case in which this could happen is one where a final distance d c (Xi, D N Pid ) 

10 between the cluster representative X b X 2 and a normalized discrete demand density D N P i d , 
D N P2d , D N P3d , D N P4d that is in the clustering step 30 determined to be a member of the 
corresponding cluster Ci, C 2 , is relatively large. A typical scenario for the clustering step 30 is 
to perform the determination of the distance dc(Xi, D N Pid ) according to a random pairing 
selection, meaning selecting from all normalized discrete demand densities D N P i d , D N P2d , D N P3d , 

15 D N P4d one after the other at random. Depending on the randomizing algorithm and also on the 
number of iterations, situations can occur in which the above described effect of a large final 
distance d c (X i5 D N Pid ) may be present. This large final distance dc(X is D N Pid ) would eventually 
lead to an undesired distance between the found continuous density D m P i c , D m P2c , D m P3c , D m P4c 
and the normalized discrete demand density D N P i d , D N P2d , D N P3d , D N P4d . 

20 A clustering correction step 60 can be used to reduce this problem. For each product Pi, P 2 , P 3 , 
P 4 a distance d p between its continuous density D m P i c , D m P2c , D m P3c , D m P4c and its normalized 
discrete demand density D N P i d , D N P2d , D N P3d , D N P4d , is determined. This is done in a threshold test 
step 62. For those of the products P b P 2 , P 3 , P 4 whose distance d p exceeds a predetermined 
threshold value t h , a substitute-selection step 63 is performed. On the condition that the 

25 threshold value t h is exceeded, the normalized discrete demand density D N P i d , D N P2d , D N P3d , D N P4d 
is enabled, depicted in form of a logical AND-function 64, as input to the substitute-selection 
step 63. 

In the substitute-selection step 63 for the corresponding normalized discrete demand density 
D N P i d , D N P2d , D N P3d , D N P4d out of the predetermined set L D of continuous model densities D m Mc a 
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substitute continuous density D m P i cS , D m P2cS , D m P3cS , D m P4c s is selected, which is defined for that 
product Pi, P 2 , P 3 , P4 as its continuous density D m P i c , D m P2c , D m P3c , D m P4c . With other words, 
instead of using the cluster-representative discrete density Rfc as input to the selection step 40, 
the selection step 40 is again performed using now as input the normalized discrete demand 
5 density D N P id, D N P2 d, D N P3d , D N P4 d itself The found normalized substitute continuous density 
D m P i cS , D m P2cS , D m P3cS , D m P4c s is used as a substitute, depicted in form of a logical OR-function 66 
in fig. 1, for the continuous density D m P i c , D m P2c , D m P3c , D m P4 c, that is only enabled, depicted in 
form of a logical AND-function 65 in fig. 1, to be used if the threshold test step 62 resulted in 
the finding that the threshold value t h is not exceeded. 

10 Once the continuous densities D m P i c , D m P2c , D m P3c , D m P4c have been determined, an adjustment 
step 70 follows, wherein for each product Pi, P 2 , P 3 , P 4 the continuous density D m P i c , D m P2c , 
D m P3c , D m P4c is adjusted into the continuous demand density D P i c , D P2c , D P3c , D P4c such that the 
continuous demand density D P i c , D P2c , D P3c , D P4c and the discrete demand densities D P i d , D P2d , 
D P3d , D P4d , respectively, have substantially the same average and the same variance. Therefor 

15 the factors s a , s v from the normalization step 20 are used, more precisely, the adjustment factor 
set Vssi, v ss2 , v ss3 , v ss4 are used. The effect is hence a reshifting and rescaling. 

Thereafter for each product Pi, P 2 , P 3 , P 4 the corresponding continuous demand density D P i c , 
D P2c , D P3c , D P4c has been determined. These continuous demand densities D P i c , D P2c , D P3c , D P4c 
can be used for determining a corresponding safety stock level S\ for each of the products Pi, P 2 , 
20 P 3 , P 4 in a safety stock determination step. In addition herefor, a given service level SL, that 
here is an ai -service level, is used. 

A possible formula for determining the safety stock level s,- is the following: 
SL(P t , s, ) = £ A (-f= + M —)p k + P 0 

with 

X 

A Pic (x)= \D Pic (t)dt 

25 
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with the lead time LT for the specific product Pi and the average |a of the continuous demand 
density D Pic of that product Pi, the probability p 0 of demand of 0 pieces of that product Pi being 
demanded during the lead time LT, and the probability p k that on a number of k days, or more 
generally time units used in the demand time series H(Pi), H(P 2 ), H(P 3 ), H(P 4 ), more than 0 
5 pieces of that product Pi are being demanded during the lead time LT. 

Generally speaking, the continuous demand density D P i C determines for the specific product Pi 
the relation between the given service level SL of that product Pj and its safety stock level Si. 
Having in addition for that specific product Pi the corresponding lead time LT, which is the 
time that passes between placing an order for for that product Pi to a product provider and the 

10 factual delivery of the product Pi to the ordering entity, the minimum quantity Q m of the product 
Pi in order to be able to guarantee the given service level SL can be determined. The minimum 
quantity Q m (Pi) is then given by the sum of the average demand for that product Pi during the 
lead time LT and the determined safety stock level Si for that product Pi. This also means that as 
soon as the current stock level S(Pi) for this product Pi falls below or reaches that minimum 

15 quantity Q m , a reordering should be initiated. This is called the reorder point. The reordering 
will then guarantee that the produces current stock level S(Pi) will not be too low. A reordering 
quantity Q r (Pi) should be typically higher than the minimum quantity Q m such that the reorder 
process and possible associated reorder costs can be held low. Since the reordering quantity 
Q r (Pi) depends on associated costs such as reordering costs, warehousing costs, depreciation, 

20 etc., it can be put at the discretion of the stock controller to determine the appropriate 
reordering quantity Q r (Pi). For the automatic reordering, therefore the system would be 
provided by the stock controller with the desired reordering quantity Q r (Pi). However, it is also 
possible to have the reordering quantity Q r (Pi) predetermined to enable an automatic reordering 
process without external intervention. 

25 There are various possible definitions of service level SL, and herein the service level SL is 
defined as the percentage of cycle times without any occurrence of a stock-out condition, i.e. 
the inability to satisfy an occurring demand for a product Pi on stock. A replenishment policy 
can be applied that triggers or initiates a product order as soon as the current stock level S(Pi) 
thereof drops to or below the minimum quantity Q m . 
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Hence, once having determined the safety stock level Si, a product order step can be performed 
wherein via an output device an order to a product Pi, P 2 , P 3 > P4 is initiated towards a 
product-providing entity, in the event that the current stock level S(Pi) of that product Pi, P 2 , P3, 
P 4 falls below the minimum quantity Q m . Thus, the method provides an automated way of 
5 controlling a product stock level under use of a determined safety stock level S\. 

The precision of the above described method can be improved by applying further exceptional 
rules. The value of 0 pieces of products Pi, P 2 , P3, P4 demanded is insofar a special value since 
it is in some sense an insignificant value. This anomaly of selective significance can be taken 
into account by the following improvements: 

10 The demand for 0 pieces of the products Pi, P 2 , P3, P4 is assigned a separate bin. 

In the normalization step 20 the demand value for 0 pieces of the products P b P 2 , P3, P4 in the 
discrete demand densities D P i d , D P2d , D P3d , D P4 d is inversely rescaled with respect to the demand 
values for more than 0 pieces of the same product Pi, P 2 , P 3 , P4. This can be done by 
determining the average zero-demand probability p a0 for the demand for 0 pieces of a product 

15 Pi, P 2 , P 3 , P 4? over all products Pi, P 2 , P 3 , P 4 , and determining the average non-zero-demand 
probability p a i+ for the demand for more than 0 pieces of a product Pi, P 2 , P 3 , P4, over all 
products Pi, P 2 , P 3 , P4. Furthermore a relation level l r is defined that gives the desired relation 
between the average non-zero-demand probability p a i+ and the average zero-demand probability 
p a0 . Then the rescaling is executed in that the demand values for 0 pieces of the products Pi, P 2 , 

20 P 3 , P 4 are each multiplied by a zero-demand rescaling factor which is the relation level l r 
divided by the average zero-demand probability p a0 , and the demand values for more than 0 
pieces of the products Pi, P 2 , P 3 , P4 are each multiplied by a non-zero-demand rescaling factor 
which is 1 minus the relation level l r divided by the average non-zero-demand probability p a i+. 

This rescaling has the effect that for slow-moving items the true demand occurrences are not 
25 outscaled and therefore reduced in their significance by the demand value for 0 pieces. 
Otherwise all slow-moving items would result to arrive in the same cluster G, although the 
density for more than 0 pieces may differ strongly. The rescaling allows to stress the significant 
values, i.e. those above 0 pieces. Although the demand value for 0 pieces is less significant, it 
should not be suppressed here, since a suppression would again distort the result of the 
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clustering. In the case of suppression, slow-moving items would be clustered together with 
non-slow-moving items, which would also not be correct, and even within the class of 
slow-moving items there is a distance between "slow" and "very slow", i.e. the products Pi, P 2 , 
P 3 , P 4 with very few demand occurrences, e.g. 3 times in a year, and those with more, e.g. 20 
5 times in a year. These distances shall not be neglected in the clustering step 30, hence the 
rescaling is seen as the appropriate measure. The rescaling hence results in a more accurate 
clustering in particular for the slow-moving items. 

In the selection step 40 the demand value for 0 pieces of the products P b P 2 , P 3 , P4 is 
suppressed. Again this is particularly useful because in comparison with the continuous model 
10 densities D m Mc the value of 0 pieces would distort the assignment of the best fitting model 
density D m Mc . Now here a total suppression has no detrimental effect, since the clustering has 
already been done. 

In the adjustment step 70 the demand value for 0 pieces of the products Pi, P 2 , P3, P4 can be 
reused and rescaled. This achieves an undistorted equivalence between the original discrete 
15 demand densities D P i d , D P2d , D P3d , D P4 d and the final continuous demand densities D P i c , D P2c , 

D P 3 C , D P 4c- 

The described method can be executed on a computer that is equipped with a computer 
program that when executed carries out the steps of the method. More precisely the computer 
would be designed to comprise a processor for executing the different steps, a memory for at 

20 least intermediately storing the results of the different steps and the predetermined set L D of 
continuous model densities D m Mc . The computer could also comprise an input device for 
receiving from a user a desired value for the service level SL and/or the reordering quantity Q r . 
The computer could furthermore be connected to a computer network and electronically effect 
the issuance of an order to a product provider for a specific product P b P 2 , P 3 , P4 at a point in 

25 time when the current stock level Si of that product Pi, P 2 , P 3 , P 4 falls below the determined 
minimum quantity Q m . The system can use a predetermined reordering quantity Q m herefor. In 
this sense the method would provide an automated solution for the problem of maintaining the 
stock of several products Pi, P 2 , P 3 , P 4 at a level that corresponds to a desired service level SL. 
The automatic order issuance can also be effected by controlling an output device of any kind, 
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e.g. printer, plotter, telephone, fax machine, loudspeaker, display, etc. whose output is 
receivable by a product provider, e.g. issuing an automatic order printout that can be sent to the 
product provider, or issuing over a phone connection a voice message that is receivable at that 
product provider, creating an order e-mail that is automatically sent to the product provider, and 
5 the like. The computer could also instead of ordering directly such a product Pi, P 2 , P 3 , P4, issue 
an alerting message to a stock controller who thereby has a stock level control system and can 
decide thereupon whether or how much of the respective product P b P 2 , P 3 , P4 is to be ordered. 

In fig. 2 a product stock control system is depicted making use of the above described method. 
The system comprises a computer 1 that receives as input for each product Pi, P 2 , P 3 , P4 its 

10 demand time series H(Pi), H(P 2 ), H(P 3 ), H(P 4 ), its lead time LT(P,), LT(P 2 ), LT(P 3 ), LT(P 4 ), its 
reorder quantity Q r (Pi), Qr(P 2 ), Qr(P 3 ), Qr(P 4 ) and the current stock level S(P,), S(P 2 ), S(P 3 ), 
S(P 4 ). It furthermore receives a service level, that can be a product-specific service level SL(Pi), 
SL(P 2 ), SL(P 3 ), SL(P 4 ). The computer executes by means of its processor the stored computer 
program that performs the safety stock determination method, determining going e.g. through 

15 the conversion step 10, normalization step 20, clustering step 30, selection step 40, parameter 
selection step 50 and denormalization step 70 for each product Pi, P 2 , P 3 , P 4 the continuous 
demand density D P i c , D P2c , D P3c , D P4c is determined. Then follows the safety stock determination 
step 80 in which for each product Pi, P 2 , P 3 , P 4 the safety stock level Si(Pi), Si(P 2 ), Si(P 3 ), Si(P 4 ) is 
determined in dependence on the product-specific service level SL(Pi), SL(P 2 ), SL(P 3 ), SL(P 4 ). 

20 In a subsequent reorder-point determination step 90, i.e. the minimum quantity Q m (Pi), Q m (P2), 
Qm(P 3 ), Q m (P4) is determined for each product Pi, P 2 , P 3 , P 4 under use of the lead times LT(Pi), 
LT(P 2 ),LT(P 3 ), LT(P 4 ). 

In a subsequent stock analysis step 100 under use of the current stock level S(Pi), S(P 2 ), S(P 3 ), 
S(P 4 ) the computer 1 determines the current stock situation for each product Pi, P 2 , P 3 , P 4 . As 
25 soon as the stock situation for at least one of the products Pi, P 2 , P 3 , P 4 becomes critical, i.e. the 
current stock level S(Pi), S(P 2 ), S(P 3 ), S(P 4 ) approaches the minimum quantity Q m (Pi), Qm(P2), 
Qm(P 3 ), Q m (P4) to a predetermined extent, e.g. becomes equal or lower than the minimum 
quantity Q m (Pi), Qm(P2), Qm(P 3 ), Qm(P4), the computer 1 will perform a stock control action for 
that product Pi, P 2 , P 3 , P 4 . The stock control action can e.g. be a stock-alert signal 101 that is 
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sent out via a computer output port. Another possibility for the stock control action is to 
perform a product order step 110 by initiating an automatic order 102 under use of the 
corresponding reorder quantity Q r (Pi), Qr(P2)> Qr(P3), Qr(P4). 

It is obvious for the person skilled in the art that the present invention can be realized in 
5 hardware, software, or a combination of these. Also, it can be implemented in a centralized 
fashion on one single computer system, or in a distributed fashion where different elements are 
spread across several interconnected computers or computer systems, whereby any kind of a 
computer system - or other apparatus adapted for carrying out the methods described herein - is 
suited. A typical combination of hardware and software could be a general purpose computer 
10 system with a computer program that, when being loaded and executed, controls the computer 
system such that it carries out the methods described herein. The present invention can also be 
embedded in a computer program product, which comprises all the features enabling the 
implementation of the methods described herein, and which - when loaded in a computer 
system - is able to carry out these methods. 

15 Computer program means or computer program in the present context mean any expression, in 
any language, code or notation, of a set of instructions intended to cause a system having an 
information processing capability to perform a particular function either directly or after either 
or both of the following a) conversion to another language, code or notation; b) reproduction in 
a different material form. 

20 Any disclosed embodiment may be combined with one or several of the other embodiments 
shown and/or described. This is also possible for one or more features of the embodiments. 

It is obvious that a person skilled in the art can modify the shown arrangements in many ways 
without departing from the gist of the invention which is encompassed by the subsequent 
claims. 



